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Using a generating-function formalism, we compute the contribution of momentum conservation 
to multiparticle correlations between the emitted particles in high-energy collisions. In particular, 
we derive a compact expression of the genuine M-particle correlation, for arbitrary M. 
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The main purpose of this paper is to deal with a rather 
general problem: consider N random variables (in a D- 
dimensional space) p 1; . . . , Pat, which are independent 
except for a constraint pi + . . . + pm — 0. What is the 
multiple correlation between M variables among the pj 
(with M < N) induced by this constraint, and in partic- 
ular, what is the cumulant of the M- variable correlation? 

One may for example think of the correlation between 
M monomers of a finite-size- N ring polymer. Another 
instance of the problem under study, from which we bor- 
row the terminology in this paper, is the constraint aris- 
ing from global momentum conservation in high-energy 
collisions: in the center-of-mass frame, the sum of the 
momenta of the N emitted particles vanishes, and this 
induces correlations between M arbitrary particles. 

An accurate knowledge of these unavoidable multipar- 
ticle correlations due to global momentum conservation is 
important. They represent an ever-present background 
effect to other sources of correlations between the out- 
going particles, which are a primary target of investiga- 
tion, as they provide information on the physics involved 
in the collisions 0: either short-range correlations, as 
quantum (Bose-Einstein or Fermi-Dirac) effects HH0, 
correlations between decay products, (mini)jets, etc.; or 
inter-particle correlations arising from a collective motion 
( "flow" ) in the context of heavy-ion collisions 5] . For in- 
stance, it has been shown that the analysis of collective 
flow can be biased by the two-particle correlations due 
to global momentum conservation Q 

A common property of short-range correlations is their 
scaling. If N denotes the total number of emitted par- 
ticles, the connected part of the M -particle correlation 
scales as 1/N l . This scaling follows from simple com- 
binatorics: 1/N 1 " 1 ^ 1 is the probability that M arbitrary 
particles have been together in a small phase-space re- 
gion. 

Global momentum conservation is certainly not a 
short-range correlation, since it affects all particles in the 
collision. One may therefore wonder how the resulting 
correlation scales with N for a given M. The two-particle 
case has already been calculated 0, Q • In the following, 
we extend the calculation to M-particle correlations for 
arbitrary M, using a generating function of the correla- 



tions and a saddle-point approximation. In particular, 
we shall derive a compact expression for the correlation 
at any order M, to leading order in 1/N, and show that 
it follows the same behavior as short-range correlations. 

We first introduce, in Sec. [I] a generating function of 
the multiparticle correlations. This generating function 
then allows us to derive the expression of the genuine 
correlation (cumulant) between an arbitrary number of 
particles to leading order in 1/N (Sec. |nj, and in par- 
ticular the scaling of the cumulants. As an illustration 
of this result, we then compute explicitly the two- and 
three-particle correlations fSec. HIIfl . The results are com- 
mented in Sec. IIVI 



I. MULTIPARTICLE CUMULANTS 

Consider a collision in which a total of N particles 
are emitted (throughout this paper, N is assumed to be 
large). Let pi, pjy, be their D-dimcnsional mo- 
menta, where the space dimension D is kept arbitrary 
for sake of generality. For a given M < N, we denote by 
f(pj 1 , ■ ■ . , Pj M ) the M-particle probability distribution 
of Pn, Pj M - _ 

To derive multiparticle correlations of arbitrary order 
in a systematic way, we introduce a generating function 
of the distributions, which we define as 

G(xi,...,x N ) = 1 + xi/(pi) + x 2 f(p 2 ) H 

+XiX 2 /(pi,P2) + • • • , (1) 
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and so on for every order M. Thus, the coefficient of the 
term x^ Xj 2 . . . Xj M is the M-particle probability distri- 
bution /(pj 1( pj a ,. • -,Pj M )- 

We shall mostly be interested in the connected part of 
the M-particle distribution, that is, the term which can- 
not be expressed as a product of correlations between less 
than M particles. For instance, the two-particle distri- 
bution can be decomposed into two parts: 

/(Pl,P2) = /(Pl)/(P2) +/c(Pl,P2), 

where the first term in the right-hand side (rhs) is the 
mere product of the one-particle probability distribu- 
tions, while the second one, the "connected part," reflects 
correlations between pi and p 2 . In particular, / c (pi,P2) 
vanishes if pi and P2 are uncorrelated. This genuine cor- 
relation, which is also called the cumulant of the distribu- 
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tion in probability theory, is obtained by taking the loga- 
rithm of the generating function of distributions • The 
cumulant fdpjt ,■■•■> Pj M ) i s the coefficient of Xj 1 . . . Xj M 
in the expansion of lnG(xi, . . . , xjv): 

lnG(^i, . . .,x N ) = x 1 f c (p 1 ) + x 2 f c (p 2 ) H 

+ x 1 x 2 f c (pi,P2) H -(2) 

If the system splits into independent subsystems or, 
more generally, if correlations in the system are short- 
range, the probability distributions of each subsystem 
add up: /({ Pj }) = J2 A (N A /N)f A ({ Pj }), where the 
sum runs over subsystems with sizes Na- It follows 
that G can be factorized into the product of functions 
for each subsystem: G({xj}) — Y\ A g a{{N ax j / N}) . In 
the limit where each individual particle is a subsystem, 
G({xj}) = [g({xj/N})] N . Therefore, InG, and the cu- 
mulants, are the sums of the corresponding quantities 
for the subsystems. As a result, inspecting the term in 
Xj t . . . Xj M shows that, while the M -particle probability 
distribution f{pj 1 , . . . , Pj M ) is independent of N, the cu- 
mulant /c(Pj!, ■ ■ • ,Pj M ) scales like l/jV M_1 . In Sec. ITT1 
we shall show that the cumulants from the correlations 
due to global momentum conservation follow the same 



scaling to leading order in 1/N. More precisely, we shall 
demonstrate that the generating function of distributions 
reads 

where the sum runs over terms with q > I. 



II. CUMULANTS IN THE LARGE- iV LIMIT 

In this Section, we derive a compact expression of the 
cumulants of multiparticle correlations due to global mo- 
mentum conservation to leading order in 1/N. In partic- 
ular, we show that the M-particle cumulants scale like 
l/N M -\ 

Let us assume for simplicity that the only source of cor- 
relation between the particles is global momentum con- 
servation, pi + • ■ • + Pat = 0. Under this assumption, the 
M-particle distribution of pi, . . . , Pm, with M < N, is 
defined as 
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' M \ „ N 

ip(Pi) s D ( P1 +-..+p N ) n [F( Pj )d D P3 ]/^ 

/(Pi,...,Pjw-) = jf , (3) 

f S D ( Pl + ■ ■ ■ + PN ) f[ [F( Pj ) d D Pj ] /ME 

where F(p) is the one-particle momentum distribution unrenormalized for the momentum conservation constraint, 
and Nd = fF(p) d D p is a normalization constant. In the following, we denote by (. . .} an F(p)-weighted average, 
that is, (g(p)) = j g(p)F(p) d D p/AfD, for any function of momentum g(p). 

The denominator in Eq. (|3J) is a constant, which we denote by 1/Cd, independent of M. The actual value will not 
influence the forthcoming discussion. 

Consider next the numerator of Eq. |J3J|. Introducing a Fourier representation of the Dirac distribution, it reads 

/ (0 (n np^--) <^>»- - / ,0 (e--)" (n ^<p^) ■ w 

Inserting in Eq. the expression of the M-particle distribution, Eq. © , with the numerator replaced by the rhs 
of Eq. (0J, we obtain 

d D k , vl V r /. , e ik 'P 



G{x u ...,x N ) = C D j (e tk 
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In passing from the first line to the second one, we have used the fact that we shall only consider the coefficient 
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of 

x ji ■ ■ - x 3m where the M indices ji, . . . , j'm are all 
different [see Eq. |J2J]. One easily checks that the only 
difference between the two forms of the generating func- 
tion in Eq. (jSJ) comes from terms in which at least one 
Xj is raised to some power m > 2 (corresponding to the 
"autocorrelation" of particle j with itself). Therefore, as 
far as we are concerned, Eq. (JSJ really is an identity, not 
an approximation. 

Please note that in Eq. (JSJ), the variable Xj is always 
multiplied by a factor F(pj). We may thus rescale Xj by 
this factor, and drop it in the following to simplify expres- 
sions. This is quite satisfactory, since it means that the 
measurable multiparticle correlations f(pj 1 , ■ ■ ■ , Pj M ) or 
f c (pj! , ■ ■ ■ , Pj M ) will not depend on the non-measurable 
distribution -F(p) — this justifies our previously calling it 
"unrenormalized" . 

To evaluate the integral in Eq. (J5J), we rely on the fact 
that N is large, and use a saddle-point approximation. 
To simplify the discussion, we introduce the notation 

^( k )=ln<e^> + £||-„. (6) 



Thus, the integrand in Eq. (0 is e N:F ( k \ We call k the 
position of the saddle point. Please note that T and its 
successive derivatives, which we shall denote by T\ T" , 
J-( 3 \ . . . , depend on x only through x/N, where x stands 
for any of the Xj. Therefore, ko, which is of course the 
solution of .F'(k) = 0, also depends on x/N only. One 
should pay attention to the fact that the saddle point ko 
is not merely the origin k = 0. 



We shall now demonstrate that to leading order in 
1/N, all multiparticle cumulants are determined by the 
saddle-point value A^JF(ko). More precisely, we show that 
the M-particle cumulant is of order 1/7V A/_1 , and that 
corrections to the saddle-point calculation only yield sub- 
leading terms, suppressed by (positive) powers of 1/N. 



Using a Taylor expansion of .F(k) around the saddle 
point, the generating function (JSJ reads 



G(x u ...,x N ) = C D e N ^«) 
= C D e N ^ 



d D k 

(2n)° 



e A^"(k )(k-k ) 2 /2 exp 
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As recalled in Sec. [I] the cumulants are given by the 
logarithm of the generating function. Now, the logarithm 
of Eq. (JJJ will split into three parts. First, there is In Co, 
which does not depend on x, and thus does not influence 
the values of the multiparticle cumulants. 

Next, there is iV.F(ko). As noted above, both T and 
ko only involve powers of x/N. Thus, the coefficient of 
Xj 1 . . -Xj M in .F(ko) contains a factor 1/N M . Multiply- 
ing by the overall factor of N, we find that the contribu- 
tion of ^^(ko) to the cumulant f c (xj 1 . . . Xj M ) scales as 
l/N M -\ 

Finally, the cumulants involve the contribution of the 
logarithm of the integral in Eq. J7J). We shall use the 
fact that the integrand is the sum of a Gaussian func- 
tion and of combinations of its moments. After integra- 
tion, this yields a common factor, l/[2irNJ 7 "(ko)] D ^ 2 , 
which multiplies a sum E. Let us show that E only 
involves terms in x l /N q with q > /. Noting that the 
values which contribute to the Gaussian integral are 
of order (k — k ) 2 ~ 1/N, one sees that each term 
A r J r(m) (k )(k - k ) m is of order N 1 -" 1 / 2 , multiplied by 
a function of x/N. Since the sum runs over m > 3, 
jyi— m/2 < ^y-i/2. therefore, a power x l goes with at 
most a factor 1/N l+1 ^ 2 . Even after raising the sum to 



the power n, and taking the logarithm, this will always 
remain of the form x l /N q with q > I. Finally, the loga- 
rithm of the overall factor l/[2-KN!F"(ko)] D ^ 2 involves x 
only through the form of powers of x/N in ln.F"(ko). All 
in all, the contribution to the M-particle cumulant of the 
integral is therefore at most of order 1/N M , subleading 
with respect to the contribution of AJF(k ). 

Therefore, we have shown that the leading contribu- 
tion to the cumulants of multiparticle correlations due 
to momentum conservation is the saddle-point value: 
In G(xi , . . . , Xji) = N!F(ka). This means first that the 
mere knowledge of k gives access to all cumulants at 
once, at least to leading order. One easily checks that ko 
is pure imaginary, i.e., iko has real- valued components. 
From Eq. it follows that .F(ko) is real- valued as well, 
and so are the cumulants, as should be. We shall illus- 
trate these points by computing explicitly the two- and 
three-particle cumulants in the following Section. Our 
result also means that genuine M-particle correlations 
arising from momentum conservation, which is a long- 
range effect, scale in the same way as correlations from 
short-range sources, namely as 0(1/A^^ /_1 ). This was 
certainly not obvious a priori. 
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III. TWO- AND THREE-PARTICLE 
CUMULANTS 

As an example of the order of magnitude derived in 
the previous Section, let us compute the two- and three- 
particle cumulants. According to the discussion in Sec. [I] 
they are given by the terms in x 2 and x 3 in the generating 
function of cumulants In G(x\ , . . . , xn), that is, following 
Sec. [H] the corresponding terms in A^ko). We shall 
for simplicity perform calculations assuming that F(p) 
only depends on the modulus |p|, not on the azimuthal 
angle of p. As a consequence, the average momentum (p) 
vanishes. Departures from this assumption are discussed 
at the end of this Section. 

Now a straightforward calculation shows that the sad- 
dle point ko is given by 



A' 



x ^ Xj < k 

2-/ 77 (e ik °-P) 



1 I (pe lko P ) 



N 

E 



^V' k "• (8) 



As mentioned above, ko is a function of x/N. Equation 
(JSJ) can be solved order by order in x/N, using the fact 
that one term (in the left-hand side) is of order in x/N 
while the other two are linear in x/N. Inspecting Eq. 10, 
and using the fact that the first term in the rhs is even 



in k, we see that the calculation of ^(ko) to order x 3 
requires our knowing ko to order x 2 (while the calculation 
to order x 2 only requires ko to order x). 

We can then compute ko , expanding Eq. (jSJ) in powers 
of ko , which gives 



ik 



A 1 



D 



D 
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D 

W) 



Xi, (9) 



where 1 d denotes the unit D x D matrix and we have 
introduced 

AT AT N 



Note that ?ko is real-valued, as expected. Reporting its 
value in Eq. ©, we obtain 



^(k ) = A - 



D 



■(Xx) 2 



2<p 2 )' 

Xi • A 1d — t^"^2 ) • Xi. 



2p 



(P 2 )" 



Multiplying this result by N, we finally obtain, to leading 
order in 1/A: 



lnC?(a;i, 



,XNj 



N 

E; 



D 



2A(p 2 



D 



•j\ E 

XH 1 3,k,l 



Pj Pi 



D 



(P 2 ) 



(Pj • Pk){pk ■ Pi 



2\ 



+ 0(x 4 ). 



Hence the coefficients of X1X2 and X1X2X3: 

Dpi ■ P2 



/c(Pl,P2) = - 
/c(Pl,P2,P3) : 



A(p2) ' 

D / X 
• jV2/p2\ ( Pl ' P2 + Pi • P3 + P2 • P3) 

D 2 

■ N 2/ p 2\2 I(Pl ' P2)(Pl ' P3j + (Pi ' P2)(P2 ' P3j + (Pi • P 3 )(P2 ■ P3)] 



(10) 



(11) 



(12) 
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We recover, in the case D = 2, the expression of the two- 
particle correlation due to transverse momentum conser- 
vation already derived in Refs 0, Q : the correlation is 
back-to-back, and stronger between high-momenta par- 
ticles. 

Finally, let us comment on the assumption that F(p) 
only depends on the modulus of p, not on its azimuthal 
angle. 

Our expansion of Eq. © relies on both (p) = and on 
the identity ((k ■ p) 2 } = k 2 (p 2 )/D, which are no longer 
valid if F is non- isotropic. In the general case, k 2 (p 2 )/£) 



will be replaced by another quadratic form in k, namely 
k • (p' ® p') • k, where p' = p — (p). To perform ex- 
plicit calculations, it becomes necessary to introduce the 
principal-axes frame, in which p' <gi p' is diagonal. For 
instance, the two-particle cumulant, Eq. will read 

'•<».»> ™ 

where the sum runs over the coordinates along the princi- 
pal axes. The distinction between the various directions 
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is for example relevant in a high-energy collision: while 
the transverse momentum distribution can be isotropic 
(unless there is some anisotropy, as, e.g., a correlation 
to the impact parameter direction), the isotropy will not 
extend to the beam direction. In that case, one may want 
to use Eq. with D = 2, when studying two-particle 
correlations in the transverse plane, but turn to Eq. 
if interested in 3-dimensional correlations. 



IV. DISCUSSION 

We have shown how it is possible to calculate the mul- 
tiparticle correlations arising from momentum conserva- 
tion between any number of particles using a generating 
function formalism: 

hxG(xi,...,x N ) =NF(k ), (14a) 

where ^(k) is given by Eq. JHJ and ko by 

•F'(ko) - 0. (14b) 

In particular, we have seen that the M-particle cumu- 
lant scales as l/iV M_1 , where N is the total number of 
emitted particles. That means that the correlation scales 
in the same way as correlations arising from short-range 
final interactions or from resonance decays, although the 
underlying reason is less obvious. In the latter cases, the 
scaling is a simple consequence of combinatorics: the M 
particles are required either to be altogether in a small 
phase-space region, or to originate from a single decay, 



hence the factor l/N^ 1 ^ 1 . However, in the case of global 
momentum conservation, the same arguments do not ap- 
ply a priori. 

The scaling of the genuine M-particle correlation due 
to momentum conservation has several consequences. 
First, a good feature: in the context of collective flow 
in heavy- ion collisions a new method of analysis has 
been proposed |9| , which relies on the idea that a cumu- 
lant expansion allows one to separate genuine collective 
phenomena from trivial short-range correlations, while 
they interfere in the data. Since correlations from mo- 
mentum conservation behave as short-range correlations, 
they can be removed as efficiently as them, to leave a 
clean collective flow signal. 

Then, some bad news. Since the M-particle correlation 
from momentum conservation is of the same order as that 
from short-range interactions, as for instance quantum 
correlations, it may bias measurements of these other cor- 
relations. It is thus worth checking that momentum con- 
servation does not contribute significantly to the correla- 
tions measured by particle interferometry and attributed 
to quantum (anti)symmetrization of the wave-function. 
This is especially true when studying the dependence of 
HBT parameters on the average momentum of the par- 
ticles since the correlation from momentum conservation 
increases with momentum [see Eqs. Ill II) and (1121) ]. 
Acknowledgments 

Fruitful discussions with Francois Gelis and Jean- Yves 
Ollitrault are warmly acknowledged. 



[1] See e.g. E. A. De Wolf, I. M. Dremin and W. Kittel, Phys. 

Rep. 270, 1 (1996) hep-ph/9508325 . 
[2] D. H. Boal, C. K. Gelbke, and B. K. Jennings, Rev. Mod. 

Phys. 62, 553 (1990). 
[3] U. W. Heinz and B. V. Jacak , Ann. Rev. Nucl. Part. 

Sci. 49, 529 (1999) nucl-th/9902020 ; U. A. Wiede- 

mann and U. W . Heinz, Phys. Rep. 319, 145 (1999) 

|nucl-th/ 9901094 ; B. Tomasik and U. A. Wiedemann, 

|nep-ph/0210250| 
[4] H. Heiselber g and A. P. Vischer, Phys. Rev. C 55, 874 

(1997) nucl-th/9609023 ; U. W. Heinz and Q. H. Zhang, 



Phys. Rev. C 56, 426 (1997) nucl-th/9701023 . 
[5] See e.g. J.-Y. Oll itrault, Nucl. Phys. A 638, 195c (1998) 

nucl-ex/9802005 . 
[6] N. Borghini, P. M. Dinh, and J.-Y. Ollitrault, Phys. Rev. 

C 62, 034902 (2000) |nucl-th/0004026]. 
[7] P. Danielewicz et al, Phys. Rev. C 38, 120 (1988). 
[8] N. G. van Kampen, Stochastic Processes in Physics and 

Chemistry (North Holland, Amsterdam, 1981). 
[9] N. Borghini, P. M. D inh, and J.-Y. Oll itrault, Phys. Rev. 

C 64, 054901 (2001) |nucl-th/0105040] . 



